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We calculate the configurational entropy S^onf for the 
SPC/E model of water for state points covering a large region 
of the (T, p) plane. We find that (i) the (T, p) dependence of 
Sconf correlates with the diffusion constant and (ii) that the 
line of maxima in Sco^f tracks the line of density maxima. 
Our simulation data indicate that the dynamics are strongly 
influenced by S'conf even above the mode-coupling tempera- 
ture rMCT(p)- 

PACS numbers: 61.43.Fs, 64.70.Pf, 65.50.-hm, 66.20.-^d 

Computer simulation studies are contributing to the 
understanding of the slow dynamical processes in simple 
and molecular liquids on approaching the glass transition 
101 . In particular, the space, time and temperature de- 
pendence of many dynamical quantities have been calcu- 
lated and compared with the predictions of the mode- 
coupling theory (MCT) which provides a descrip- 
tion of the dynamics in weakly supercooled states [|| of 
many atomic and molecular liquids, including the SPC/E 
model of water studied here |Qj5| . 

Simulations have also been used to investigate the 
"thermodynamic approach" to the glass transition , 
which envisages a relation between diffusion constant D 
and configurational entropy S'conf by relating the dynam- 
ics of liquids at low T to the system's exploration of 
its configuration space. The properties of the liquid are 
dominated by the basins of attraction of local potential 
energy minima; the liquid experiences vibrations local- 
ized around a basin and rearranges via relatively infre- 
quent inter-basin jumps Motivated by the separa- 
tion of time scales, the total entropy S may be sepa- 
rated into two parts: (i) an intra-basin contribution S'vib, 
which measures the vibrational entropy of a system con- 
strained to reside within a basin, and (ii) a configura- 
tional contribution, S'conf, quantifying the multiplicity of 
basins [ p^ . Explicit calculation of Sconf has been per- 
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formed for hard spheres soft spheres 0, Lennard- 
Jones systems and tetravalent network glasses |l^] 

with the aim of evaluating the Kauzmann temperature 
Tk, at which S'conf appears to vanish p^ . 

One of the most studied models of molecular liquids is 
the SPC/E potential, designed to mimic the behavior of 
water |15|] . The dynamical properties of this model have 
been studied in detail in the weakly supercooled regime 
and have been shown to be consistent with the predic- 
tions of MCT The SPC/E potential is of particu- 
lar interest for testing theories of the supercooled-liquid 
dynamics, since, as observed experimentally for water, 
along isotherms D has a maximum as function of the 
pressure P or of the density p [|l6|. This maximum be- 
comes more pronounced upon cooling. Furthermore, the 
line of isobaric density maxima, which corresponds to 
the line of isothermal entropy maxima from the Maxwell 
relation {dV/dT)p = —{dS/dP)T, appears strongly cor- 
related to the line of D maxima [Fig. |l|(a)], providing 
motivation to test possible relationships between Sconf 
and D. A test of such relationship, based on the analysis 
of experimental data at ambient pressure, was performed 
by Angell and coworkers in 1976 [ pTf . 

Here we calculate the entropies S, S'vib and S'conf for 
the SPC/E potential for state points covering a large 
region of the (T, p) phase diagram. We then compare 
the behavior of Sconf with D. We also estimate Tk and 
consider its relation with the MCT temperature Tmct- 
We first calculate S for a reference point (T — 1000 K 
and p = 1.0 g/cm^) ||l8|, and then calculate the entropy 
as a function of (T, p) via thermodynamic integration |ig| ] 
using the state points simulated in Ref. (s) , as well as new 
simulations extending to higher temperatures We 
show an example of S along the p = 1.0 g/cm'^ isochore 
in Fig.|(a). 

We calculate Svib, i-e. the entropy of the liquid con- 
strained in a typical basin, using the properties of the 
basins visited in equilibrium by the liquid. For each 
(T, p), we calculate, using the conjugate gradient algo- 
rithm [2l| ], the corresponding local minima - called in- 
herent structures (IS) - for 100 configurations [ p3[ . 
We estimate S'vib by adding anharmonic corrections to 
the harmonic contribution S'harm p^ ] 
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where Wi are the normal-mode frequencies of the IS de- 
termined form the Hessian matrix Ref. [|l2| found 
that the harmonic approximation for a binary Lennard- 
Jones mixture is a vahd estimate of 5vib for temperatures 
around Tmct- However, in the case of the SPC/E po- 
tential, we find that there are significant anharmonicities 
in the basins. This can be seen from the fact that, if 
the system were purely harmonic, the energy E should 
equal the energy of a minimum i^conf plus the contri- 
bution -Eharm — (6A^ — 3)/cbT/2 of the harmonic solid 
approximation. In contrast with the binary mixture 
Lennard-Joncs case, we find the vibrational contribution 
i?vib = E — i?conf to be roughly 10% larger than the har- 
monic approximation, even at the lowest temperatures 
studied. We then estimate the anharmonic contributions 
to S'vib by heating the IS at constant volume and mea- 
suring the deviation of -Evib from the harmonic approxi- 
mation [^6| . For each collection of basins corresponding 
to a particular (T, p) state point, we calculate -Evib in 
the range T = — 200 K and fit our results using the 
approximation 
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By integrating the relation dS^ib — dE^ih/T, wc find 
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S'vib = S'harm + SflT -|- -6T^ (3) 

We use Eq. (|^) to extrapolate Svib to higher tempera- 
tures. 

^From the knowledge of Svib, we calculate Sconf = S — 
Svib- As an example of the entire procedure, we show in 
Figs. 11(a) and H(b) the T dependence of Svib and Sconf 
for the yO = 1.0 g/cm'^ isochore. In Fig. ^(a), we show 
Sconf as a function of p for several T, and in Fig. ||(b) 
we show the behavior of D along the same isotherms. 
Fig. ^ shows that at high densities the behavior of ln{D) 
is nearly linear when plotted as a function of (TSconf)"^, 
as proposed by Adam and Gibbs ^ over 30 years ago. 

Fig. 0(b) shows the lines of maxima for D and Sconf in 
the region T < 260 K, where we have been able to clearly 
detect a maximum. To highlight the differences between 
S and Sconf, we show also the line at which S{p, T) has 
a maximum. Fig. ^(b) shows that the lines of maxima in 
D and in Sconf track each other within the uncertainties 
of our calculations. 

Note that the density where D has a maximum (p « 
1.15 g/cm"^), as well as the density where the maximum 
in Sconf occurs, depends weakly on T. One possible ex- 
planation is that, for the values of p where the maxima in 
D and Sconf occur, two density-dependent (and primarily 
T independent) mechanisms balance. Increasing p in wa- 
ter from the "ideal" tetrahedral density (i.e., that of ice 



Ih, about 0.92 g/cm^) leads to the progressive destruc- 
tion of the hydrogen-bond network, and hence increases 
the number of minima (and therefore Sconf) , since there 
are more configurations corresponding to a disordered 
tetrahedral network. However, at large enough density, 
core repulsion begins to dominate the liquid properties, 
as expected for "typical" liquids. In such a case, one ex- 
pects that increasing p decreases the number of minima 
in the potential energy landscape, as the system becomes 
more densely packed, and hence fewer configurations are 
possible. 

We observe that the close connection between D and 
Sconf shown in Fig. ^ occurs in the same region where the 
dynamics of SPC/E water can be rather well described by 
MCT, suggesting that MCT may be able to capture the 
reduction of the mobility due to entropic effects. More- 
over, in the same region of the (T, p) plane, D also corre- 
lates well with the number of directions in configuration 
space connecting different basins |^^. This suggests the 
possibility of a statistical relation between the number of 
minima and their connectivity ||2^ . 

Before concluding, we note that the present approach 
allows us to estimate the locus of the Kauzmann temper- 
ature Tk{p) at which Sconf would disappear on extrap- 
olating to lower T. Fig. || shows the calculated values 
of Tk{p) together with the locus of mode coupling tem- 
peratures Tmct{p)- The ratio Tmct/Tk has been used 
as an indication of the fragility of a liquid [Q; we find 
Tmct/Tk ~ 1.05-1.15 suggesting that SPC/E is an ex- 
tremely fragile liquid in the region of temperatures near 
Tmct — as also found in experimental measurements. 
We note that the values of Tk (p) depend strongly on the 
validity of the extrapolation of the potential energy pC| ] 
to temperatures lower than the one we can equilibrate 
with the present computer facilities. Slower changes in 
i?conf below the lowest simulated state point would pro- 
duce a much slower decrease for Sconf |3l[| . This may be 
quite plausible since it appears that i?conf is approach- 
ing the crystal value, which is expected to always be less 
than the liquid value. A much slower change in Sconf 
is expected to be accompanied by a slower decrease of 
D, which may be related to a possible fragile-to-strong 
transition in water ]3^ . 
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FIG. 1. (a) Line of maxima of the entropy 5* and of the dif- 
fusion constant D from experimental data on water [Q. (b) 
Lines of entropy maxima S, configurational entropy maxima 
Sconf a-nd diffusion constant maxima D (from Ref. |5[) for the 
SPC/E model. 
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FIG. 2. Calculated values for the p — 1.0 g/cm'^ isochore 
of: (a) S determined from simulation (filled □). The line is a 
fit used to extrapolate S to lower T [B0|. (b)S'yii5. (c) S'confi 
obtained by methods described in the text. 
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FIG. 3. (a) The configurational entropy Sconf along 
isothermal paths; from top to bottom T — 300, T — 260, 
T = 240, T = 230, T = 220, T = 210 K. (b) Diffusion 
constant D along the same isotherms. Notice the close corre- 
spondence of the maxima of Sconf and D. 
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FIG. 4. Diffusion constant D versus (TSconf) ^ for all 
studied isochores. 
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